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Abstract 

The completely bounded trace and spectral norms in finite dimensions are shown to be 
expressible by semidefinite programs. This provides an efficient method by which these norms 
may be both calculated and verified, and gives alternate proofs of some known facts about 
them. 

1 Introduction 

Linear mappings from one space of operators to another, which are often called super-operators, 
play an important role in quantum information theory. Quantum channels in particular, which 
model general discrete-time changes in quantum systems, are represented by super-operators that 
act on operators on finite-dimensional complex vector spaces. 

It is natural to consider distances between quantum channels, so as to quantify the similarity 
with which they act on quantum states. One way to define such a notion is to define a suitable 
norm on the space of super-operators in which channels of a given size are represented. Then, 
the distance between two channels is defined as the norm of their difference. A natural question 
that arises is: what norms give rise to the most physically meaningful notions of distance? As is argued 
in BGLN05) , the answer to this question may depend on the problem at hand — but perhaps the 
most natural and widely applicable choice within quantum information theory is the completely 
bounded trace norm, also known as the diamond norm. This norm was first used in the setting of 
quantum information by Kitaev [Kit97j, who used it mainly as a tool in studying quantum error 
correction and fault-tolerance. It is equivalent, up to taking the adjoint of a super-operator, to 
its spectral norm variant, which is usually known simply as the completely bounded norm. The 
completely bounded norm, as well as variants that include the completely bounded trace norm, 
have been studied in operator theory for many years. (See |Pau02J for historical comments and 
further details.) 

The definition of these completely bounded norms, which can be found in the section follow- 
ing this introduction, may seem unusual at first glance. It turns out, however, that they are quite 
natural and satisfy many remarkable properties. They are, in particular, much easier to reason 
about and to work with than the seemingly simpler super-operator norms that are induced by 
the trace norm and spectral norm, primarily because the completely bounded norms respect the 
structure of tensor products while the induced norms do not. The physical importance of this 
property, within the setting of quantum information theory, has been discussed in several sources 
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tional references that highlight the properties and uses of completely bounded norms in quantum 
information include |D]'KR06t |jen06| IPGWP+08I . 



One obvious question that comes to mind about the completely bounded trace and spectral 
norms is: can they be efficiently computed? Unlike the norms of operators that are most typically 
encountered in quantum information theory which are trivially computable from spectral or 
singular-value decompositions, the computation of completely bounded norms is not known to 
be straightforward. To the author's knowledge there are only two papers written prior to this one, 
namely !Zar06| and [JKP09 1, that present methods to compute the completely bounded spectral 



or trace norm of a given super-operator. Both papers describe iterative methods, and analyze the 
complexity of each iteration of these methods, but do not analyze their rates of convergence. So, 
although these papers may describe potentially efficient methods, they do not include complete 
proofs of their efficiency. 

The purpose of this paper is to explain how the completely bounded trace norm of a given 
super-operator (and therefore its completely bounded spectral norm as well) can be expressed as 
the optimal value of a semidefinite program whose size is polynomial in the dimension of the 
spaces on which the super-operator acts. Using known polynomial-time algorithms for solving 
semidefinite programs, one obtains a provably efficient algorithm (meaning a deterministic poly- 
nomial time algorithm) for calculating these norms. This approach also has the obvious practical 
advantage that it is more easily implemented through the use of existing semidefinite program- 
ming optimization libraries, and allows one to take advantage of the extensive work that has been 
done to solve semidefinite programs efficiently and accurately. Moreover, through semidefinite 
programming duality, one obtains a means by which a certificate of the value of the completely 
bounded trace or spectral norm of a given super-operator may be quickly verified. 

In a recent paper written independently from this one, Ben-Aroya and Ta-Shma [BATS09J have 
found a different (but related) way to efficiently compute the completely bounded trace norm of 
super-operators using convex programming. 

The essence of the semidefinite programming formulation of the completely bounded trace 
norm that is described in this paper appears, at least to some extent, in the paper [KWOO]; al- 
though it was not made explicit or considered in full generality therein. The present paper aims 
to present this formulation explicitly and without any discussion of the quantum interactive proof 
system model of computation, which is the primary focus of [KWOO J. A second semidefinite pro- 
gramming formulation of the completely bounded trace norm is also presented, based on the 
competitive quantum game framework of [GW07J. This formulation is slightly simpler, but is valid 
only for super-operators that are the difference between two quantum channels — which happens 
to be an important special case in quantum information. 

Semidefinite programming is useful not only as a computational tool, but as an analytic tool as 
well. The last section of this paper gives two examples along these lines that are derived from the 
more general semidefinite programming formulation of the completely bounded trace norm. The 
first example concerns an alternate characterization of the completely bounded trace norm and 
the second illustrates a precise sense in which two known characterizations of the fidelity function 
(given by Uhlmann's Theorem and Alberti's Theorem) are dual statements to one another. 

2 Background 

The two subsections that follow aim to provide the reader with an account of the background 
knowledge assumed in the remainder of the paper. The first subsection discusses well-known 
concepts from finite-dimensional operator or matrix theory, and is mainly intended to make clear 
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the notation and terminology that is used later. It also includes the definitions of the completely 
bounded norms that are the focus of this paper. The second subsection discusses semidefinite 
programming, using a form that is less common than the so-called standard form of semidefinite 
programs, but that is equivalent and better suited to the needs of this paper. 

Operators and super-operators 

The scripted letters X, y , Z, and W will denote vector spaces of the form C" for n > 1, whose 
elements are identified with n-dimensional column vectors. The j-th elementary unit vector in 
such a space is denoted e,-, and the unit sphere in X (with respect to the Euclidean norm) is denoted 

S(X) = {u G X : \\u\\ = 1}. 

For X = C" and y = C m , the space consisting of all linear operators of the form A : X — > y 
is denoted L (X, y) and is identified with the set of m x n complex matrices in the usual way. 
The notation L (X) is shorthand for L (X, X), and the space X is identified with L (C, X) when 
necessary. An inner product on L {X , y) is defined as {A, B) = Tr(A*£>) for all A, B G L (X , y), 
where A* G L (y, X) denotes the adjoint (or conjugate transpose) of A. The identity operator on 
X is denoted lx, and for each choice of indices i,j we write Ey = 

Three operator norms are discussed in this paper: the trace norm, Frobenius norm, and spectral 
norm, defined as 

II AIL = TrVA*A, \\A\L = J(A,A) , and II AIL = max ||Aw|| , 

111 V MCO ue s(X) 

respectively, for each A G L (X,y). For every operator A it holds that 

||AL < ||A|| 2 < \\A\\ X . 
The following special types of operators are also discussed. 

1. An operator X G L (X) is Hermitian if X = X* . The set of such operators is denoted Herm (X). 

2. An operator P G L (X) is positive semidefinite if it is Hermitian and all of its eigenvalues are 
nonnegative. The set of such operators is denoted Pos (X). The notation P > also indicates 
that P is positive semidefinite, and more generally the notations X < Y and Y > X indicate 
that Y — X > for Hermitian operators X and Y. 

3. An operator P G L (X) is positive definite if it is both positive semidefinite and invertible. The 
set of such operators is denoted Pd (X). The notation P > also indicates that P is positive 
definite, and the notations X < Y and Y > X indicate that Y — X > for Hermitian operators 
X and Y. 

4. An operator p G L (X) is a density operator if it is both positive semidefinite and has trace equal 
to 1. The set of such operators is denoted D (X). 

5. An operator U G L (X) is unitary if U*U = lx- The set of such operators is denoted U (X). 

A super-operator is a linear mapping of the form O : L (X) — > L (y), and the space of all 
mappings of this form is denoted T [X , y). The identity super-operator on L [X) is denoted \(x)- 
The adjoint super-operator to O G T (X, y) is the unique super-operator O* G T ( y, X) for which 
(Y,0(X)) = (0*(Y),X) for all X G L (X) and Y G L (y). The following special types of super- 
operators are discussed. 
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1. <£> G T (X, y) is Hermiticity-preserving if O(X) G Herm ( J>) for every X G Herm (A"). 

2. $6T (A", is completely positive if it holds that 

(<t>®l L(W) )(P) G Pos(3^0 W) 

for every choice ofW = C k and P G Pos (AT ® W). 

3. O G T (X, y) is trace-preserving if Tr(0(X)) = Tr(X) for every X G L (X). 

4. O G T ( X, y) is a quantum channel if it is both completely positive and trace-preserving. 

The Choi-Jamiotkowski representation J(<£>) G L (y ® A") of a super-operator O 6 T (A", J 7 ) is the 
operator defined as 

/(0)= £ «>(Ei,;)®Ei,/ 
l<i,j<n 

(where this expression assumes A" = C"). The mapping / is a linear bijection from T (X,y) to 
h(y ® X). The operator /(O), written as an nm x hih matrix, represents one convenient way that 
a super-operator may be expressed in concrete terms. 

A pair of operators (A, B) in L (X, y ® Z) is a Stinespring pair for O G T (A", J 7 ) if it holds that 

<D(X) =Tr 2 (AXB*) (1) 

for all X G L (X), and an expression of the form Q]) is called a Stinespring representation of 
Stinespring representations exist for every super-operator, provided the space Z has dimension at 
least rank(/(0)). It is straightforward to compute such a Stinespring pair (A, B) from the Choi- 
Jamiolkowski representation of O G T (X, y): for any expression 
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it holds that 

i,j,i 

forms a Stinespring pair of O. 

For every super-operator O G T (X, y), one defines the induced super-operator norms: 

\\<$>\\ x =max{||0(X)|| 1 : XeL(X), {{X^ < 1}, 

||O|| co = max{||O(X)|| 0O : X G L (X) , ||XL < 1}, 

as well as completely bounded variants of these norms: 



sup 

k>\ 



O® 1 



L(C*) 



and 



sup 

k>l 



<5®1 



L(C*) 



As was done in the introduction, we will refer to ||| <3> ||| a as the completely bounded trace norm and to 
lll^lll oo as the completely bounded spectral norm. It is common that ||| < ^ > ||| a is denoted || <J> || and called 
the diamond norm, and that ||| *!* ||| oo is denoted || O || cb and called simply the completely bounded norm. 
It holds that 



o®i 



and 



O® 1 



and that 1 1 1 <3> 1 1 1 a 
respect to tensor products, meaning that 



I ^ III oo ^ or every O G T(X,y). These norms are both multiplicative with 



l<t>®T|| 



lOli 



|T| 



and 



IO® Til 



1*1' 



ITI 



for any choice of super-operators O and T. 
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Semidefinite programming 

This section gives a brief overview of semidefinite programming, which is discussed in greater 
detail in several sources (including |Ali95. VB96, Lov03, dK02], for instance). The particular for- 
mulation that is described here is somewhat different than the well-known standard form that is 
used by most authors, but it is equivalent and more convenient for the purposes of this paper. 
A semidefinite program over X = C" and y G C m is specified by a triple (Y, A, £>), where 

1. Y G T (X, y) is a Hermiticity preserving super-operator, and 

2. A G Herm (X) and B G Herm (y) are Hermitian operators. 

The following two optimization problems are associated with such a semidefinite program: 
Primal problem Dual problem 

maximize: (A,X) minimize: (B,Y) 

subject to: Y(X) < £>, subject to: Y*(Y) > A, 

XePos{X). YGPos(^). 

With these problems in mind, one defines the primal feasible set A and the dual feasible set B as 

A = {X G Pos (X) : Y(X) < £>} , 
B = {Y G Pos (y) : Y*(Y) > A} . 

Operators X G A and Y G B are also said to be primal feasible and dual feasible, respectively. For the 
sake of the discussion of computational efficiency below, it will be helpful to also define, for each 
£ > 0, the sets 

A e = {X G Pos (X) : X + H G A for all H G Herm (X) satisfying ||H|| 2 < e} , 
B E = {Y G Pos (y) : Y + HGB for all H G Herm (y) satisfying ||H|| 2 < e} . 

Intuitively speaking, A e contains primal feasible operators that are not too close to the boundary 
of the primal feasible set, and likewise for B e . 

The functions X i— > (A, X) and Y i— > (B,Y) are called the primal and dual objective functions, 
and the optimal values associated with the primal and dual problems are defined as follows: 

a = sup {A, X) and & = inf (B, Y) . 

(If it is the case that A = or B = 0, the above definitions are to be interpreted as a = — oo and 
B = oo, respectively.) The supremum and infimum cannot always be replaced by the maximum 
and minimum — in some cases even finite values a. and B may not be achieved for any choice of 
X G A and Y G B. 

Semidefinite programs have associated with them a powerful theory of duality, which refers 
to the special relationship between the primal and dual problems. The property of weak duality, 
which holds for all semidefinite programs, is stated in the following theorem. 

Theorem 1 (Weak duality). For every semidefinite program (Y, A, B) as defined above, it holds that 
a < B. 
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This property implies that every dual feasible operator Y G B provides an upper bound of (£>, Y) 
on the value {A, X) that is achievable over all choices of a primal feasible X G A, and likewise 
every primal feasible operator X G A provides a lower bound of (A, X) on the value (B, Y) that is 
achievable over all choices of a dual feasible Y G B. 

It is not always the case that oc = B for a given semidefinite program (Y, A, £>), even when <x 
and B are finite. For most semidefinite programs that arise in practice, however, it is the case that 
a = B, which is a situation known as strong duality. There are different conditions under which 
this property is guaranteed, one of which is given by the following theorem. 

Theorem 2 (Slater-type condition for strong duality). The following two implications hold for every 
semidefinite program (Y, A,B) as defined above. 

1. Strict primal feasibility: If B is finite and there exists an operator X > such that Y(X) < B, then 
a. = B and there exists Y G B such that (B, Y) = B. 

2. Strict dual feasibility: If a. is finite and there exists an operator Y > such that Y*(Y) > A, then 
cc = B and there exists X G A such that (A, X) = a. 

One may consider a general computational problem that asks for the optimal primal and dual 
values of a given semidefinite program, possibly up to some specified accuracy. There are various 
ways in which this may be done, one of which is to phrase the problem as a promise problem [ESY84J 
such as the following one. 

Problem 3. The semidefinite programming approximation problem is as follows. 

Input: A semidefinite program (Y, A, B) over X = C and y = C m , an accuracy parame- 
ter e > 0, and a positive integer R. 

Promise: The set A £ is non-empty, and for every X G A it holds that ||X|| 2 < R- (In the 
terminology of flGLS93fl , the primal feasible region A of (Y, A, B) is well-bounded, 
with parameters e and R.) 

Output: A real number 7 such that 1 7 — a | < e, where a is the optimal value of the primal 
problem associated with (Y, A, B). 

The description of this problem does not explicitly state how the super-operator Y is to be repre- 
sented, but we will assume it is specified by the matrix representation of /(Y). Other forms, in- 
cluding Stinespring representations and Kraus representations, are easily converted to this form. 
It is also assumed that the entries of /(Y), A, and B have rational real and imaginary parts. 

The computational problem stated above can be solved in polynomial time using the ellipsoid 
method JGLS93), as the following theorem states. 

Theorem 4. There exists an algorithm that solves the semidefinite programming approximation problem 
stated above that runs in time polynomial in n, m, log(R), log(l/e), and the maximum bit-length of the 
entries of /(Y), A, and B. 

Here, the bit length of a complex number z = (a/b) + i(c/ d) is the number of bits needed to 
represent the 4-tuple (a, b, c, d), where a, b, c, and d are integers represented in binary. 

Note that the above problem asks only for an approximation to the optimal primal value, but 
the simple transformation (Y, A, B) — > (— Y, — B, —A) shows that any algorithm for it also allows 
one to approximate the optimal dual value. (Alternately, the ellipsoid method can be applied 
directly to the dual problem.) 
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It is possible to approximate more general classes of semidefinite programs efficiently. For in- 
stance, the bound || X|| 2 < R need not hold for every primal feasible X, provided certain assump- 
tions are known about the size of the optimal solution. These generalizations are not important 
for this paper, and the above problem can be more easily fit to the general presentation of [GLS93J 
(which is described in the setting of semidefinite programming in [Lov03j). 

It should be noted that one would typically not use the ellipsoid method to solve semidefinite 
programming problems in practise, given that interior point methods [Ali95, dK02J are significantly 
faster. In strictly formal terms, however, interior point methods have not been proved to run in 
polynomial time using the Turing machine model of computation. 

3 A semidefinite program for the completely bounded trace norm 

We will now describe and analyze a semidefinite program whose optimal (primal and dual) value 
is III^HIp where OeT (X , y) is an arbitrary super-operator given by a Stinespring representation 

<E>(X) = Tr z (AXB*) 

for A, B G L [X, y (g) Z). It is assumed further that Z has the minimal dimension dim(Z) = 
rank(/(0)) for which such a Stinespring representation exists. 

The primal and dual problems for the semidefinite program we will consider may be stated 
informally as follows: 

Primal problem Dual problem 

maximize: (BB*,W) minimize: || A*(ly <g> Z)A|| 00 

subject to: Tr y (W) = Try (ApA*) , subject to: ly <g) Z > BB* , 

p&D(X), ZgPos(Z). 

w g Pos (y®z). 

These problems are associated with the semidefinite program that is more formally specified as 
follows. We define a Hermiticity-preserving super-operator 



as 

Y 

The adjoint super-operator 
is given by 



y : L(xe {y®z)) -> L(cez) 

X • \ _ /Tr(X) 
• WJ ~ \ Try(W — AX A*) 

Y* : L (C © Z) -» L [X © (y ® Z)) 
A -\ _ (Kl x - A*(t y ®Z)A 



^ ' ZJ \ ty<S)Z 

(In these expressions of Y and Y*, the symbol • denotes an operator or vector of the appropriate 
dimensions upon which the output of these super-operators does not depend.) We also define 
C G Herm (X ffi (y ® Z)) and D G Herm (C © Z) as 

C= (o BB*) and D= (l 
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Now, the primal and dual problem associated with (Y, C, D) may be expressed as follows: 



Primal problem 



maximize: (BB*, W) 
subject to: Tr y (W) < Try (AX A*) , 
Tr(X) < 1, 

x e Pos(^), 
w g Pos (y®z). 

Notice that for any choice of a primal feasible operator 



Dual problem 
minimize: A 
subject to: Xl x > A*(ly 
ty ® Z > BB*, 
A > 0, 

Z G Pos (Z) . 



Z)A 



X 
M* 



M 
W 



(2) 



there exist operators P G Pos (X) and Q G Pos O <g> Z) such that Tr(X + P) = 1 and 

Tr y (W + Q) = Try (A(X + P)A*) . 

The operator 

fX + P M \ 
\ M* W + QJ 

is therefore primal feasible, and obtains at least the value achieved by (O (by virtue of the fact 
that BB* is positive semidefinite). This accounts for the informal statement of the primal problem 
above, where the inequality constraints are replaced by equality constraints. The dual problem 
above is obviously equivalent to its informal statement, because A* (ty (g) Z)A is positive semidef- 
inite for positive semidefinite Z, and therefore 



min {A > : Alx > A*(ty ® Z)A} = \\A* (ty ®Z)A\ 



Strong duality 

We will first verify that strong duality holds for the above semidefinite program, using Theorem[2] 
First note that it is clear that the optimal primal value a. is finite, for it must hold that Tr(W) < 
II ^ IlL f° r an y prirrml feasible operator with the form and therefore a < 11^)1^ || B ||^. 

Now, to verify strict dual feasibility, suppose that }i and A are positive real numbers such that 
]i > || B |£ and A > pi \\ A |£. Then 

(X \ # fX \ _ (Xt x - ]iA*A \ fO \ 

[O H t z ) >{} and ^ [o j0. z )-\ } tty®t z ) > {0 BB*)' 

which illustrates strict dual feasibility. Thus, by Theorem |2l the optimal value a associated with 
the primal problem is equal to the optimal dual value B, and is achieved for some choice of a 
primal feasible operator. 

One may wonder whether the semidefinite program above is also strictly primal feasible. Hav- 
ing already established strong duality, it is not really essential that this is proved, but it may 
be noted that strict primal feasibility indeed does hold, relying on the assumption dim(Z) = 
rank(/(0)). This observation, which happens to imply that the optimal dual value is achieved for 
some dual feasible operator, will follow from the discussion of computational efficiency below. 



8 



Optimal value 

Now let us verify that the optimal value oc = B of our semidefinite program is equal to 1 1 1 <I> 1 1| ^ - 
Define W = C k for k = max{dim(A') / dim(^ ®Z)}. Given that dim(W) > dim(AT), it holds that 

|||3>|||? = max \(U,Tr z ((A®l w )uv* (B* tg>l w )))| 2 

u,v<=S(X®W) 
UeU(V(8)W) 

= max \\(B* ®1 W )(LT* ®t z )(A®t w )u\\ 2 

ueS(x®W) 

UGU(Vcx)W) 

= max u*(A*®l w )(U®l z )(BB*®lw)(W®l z )(A®l w )u 

u<=S(X®W) 

UG\j(y(g!W) 

= max ( BB* ,Tr w [(U* ®t z ) (A ®t w )uu*( A* ®t w )(U®t z )}). 

ueS{X®W) 
UGU(Vcx)W) 

Now define two sets Q, TZ C Pos (3^ ® -2) as 

Q = {WG Pos (J ® Z) : Trj;(W) = Txy{ApA*) for some choice of p G D (A")} , 

ft = {Tr w [(Ji*(g)l 2 )(A«)l >v ) MM *(A*«)l >v )(!J(g)l 2 )] : u G <S(A' W) , U G U (3> ® W)} . 

Our interest in the set TZ is clear, for the equation above has established that 

= max (BB* , W) . 

The set Q, on the other hand, is of interest because the optimal value a. of the primal problem for 
the semidefinite program defined above is given by 

a = max (BB*, W) . 

To establish that a. = if therefore suffices to prove that Q = TZ, which is easily done as 

follows. 

First consider an arbitrary choice of u G S(X ® W) and U G U (y ® W), and let 

W = Tr w [(IT ® 1 2 )(A <g> 1 W )« M *(A* ® 1 W )(IZ ® l z )] . 

Then Try (W) = Try (A Tr w (uu*)A*), and so it holds that W G Q, which proves KCQ. 

Now consider an arbitrary element We Q, and let p G D (X) be a density operator satisfying 
Try(W) = Try (ApA*). Given that we have chosen W to have dimension at least as large as that 
of both X and y ® Z, there must exist vectors u G S(X ® W) and a? G J'^^^W such that 
p = Tryy(«M*) and W = Tryy(ow*). This implies that 

Try^w(zvw*) = Try rS)W ((A ® l w )uu* (A* ® lyv)), 

so there must exist U G U (y ® W) such that (IP <g> lg) (A <g> l w )u = w. Therefore 

W = 7r w (ww*) = Tr w [(LZ* <g> 1 2 )(A <g> l w )wu*(A* ® l w ) (U ® lg)] , 

which proves that W £ 1Z, so that Q C Kas required. 
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Computational efficiency 

Now let us verify that the optimal value W^HIi of the semidefinite program described above can 
be approximated by an efficient computation. By Theorem H] our task is to argue that suitable 
parameters R and e for the promise in Problem|3]can be determined. 

For the sake of clarity, let us summarize our notation: we have X = C n , y = C m , and Z = C , 
and OgT ( X, 3^) is the super-operator given by 

<E>(X) = Tr z (AXB*) 

for which we wish to approximate IH^HIi- The semidefinite program that represents this quantity 
is represented by the Hermiticity-preserving super-operator Y G T (X © (y © Z) ,C © Z) and 
Hermitian operators C G Herm (X © (y © Z)) and D G Herm (C © Z) as described above. We 
will take N to be the total bit-length of this semidefinite program, which is polynomially related 
to n, m and the maximum bit-length of the entries of A and B. 

First, it is clear that every primal feasible operator has trace bounded by 1 + || A \\ 0O . Given that 
the Frobenius norm is upper-bounded by the trace for positive semidefinite operators, it therefore 
suffices to choose R = 1 + || A ||^, which is obviously bounded by 2 cN for some positive integer 
constant c. 

The specification of e is slightly more complicated. Consider first the operator Try (A A*). We 
have chosen Z to have minimal dimension to admit a Stinespring representation of O, and from 
this assumption it follows that Try (A A*) is positive definite. Using the assumption that the real 
and imaginary parts of the entries of A are rational, along with the fact that nonzero roots of 



integer polynomials cannot be too close to zero (see, for instance, Theorem 2.9 of |Bug04|), one 
may derive a lower-bound on the smallest eigenvalue of Try (A A*). For the purposes of this 
analysis, it suffices to note that there exists an integer constant do > 1 such that for 5 = 2~ d ° N we 
have that the smallest eigenvalue of Try (AA*) is at least S, and therefore Stz < Try(AA*). 
Now consider the operator 

'X 
,0 W, 

where 

X = ^-l x and W= -^ty®Try(AA*), 
4n 8nm y yy ' 

along with any choice of a real number e > that satisfies 

e < 



8nm 

Let us note, in particular, that this bound holds for e = 2~ dN for some choice of a positive integer 
constant d. It is our goal to show that every Hermitian operator whose distance from P is at most 
e (with respect to the Frobenius norm) lies within the primal feasible set A, and therefore that A £ 
is nonempty. In other words, for any choice of operators H G Herm (X), K G Herm (y © Z), and 
M G L (y © Z, X) satisfying 

H M N 
M* K 



< e, 

2 



we wish to prove that 

'X + H M 



M* W + K' (3) 
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is primal feasible. 

It is clear that el < P, and therefore (O is positive semidefinite. As || K < e it follows that 

W + K < W + ety 0Z < —l y ®Tr y (AA*) 

and therefore 

7ry(W + K)<±-Jr y (AA*). 

As IIHH^ < e it holds that 

-U* < X - el x < X + H 
2n 

and therefore 

±-Tr y {AA*) <Tx y (A(X + H)A*). 

It follows that Tr;y(W + K) < Tr y (A(X + H)A*) and therefore the above operator (0 is primal 
feasible as required. 

We have shown that the requirements of the promise in Problem [3] are met for R = 2 cN and 
e = 2~ dN for some positive integer constants c and d. By Theorem|4]the value 1 1 1 <I> 1 1| may therefore 
be approximated to within error e in time polynomial in n, m and the size of the entries of A and 
B. (It is possible of course to choose a smaller error, e = 2~p( n ^ for any polynomial p for instance, 
if this is desired.) 



4 A simpler semidefinite program for quantum channel distance 



A somewhat simpler semidefinite program exists for the completely bounded trace norm of the 
difference between two quantum channels, which is a special case that is relevant to quantum 
information. This case was discussed in [GLN05J, and shown to reduce to a convex optimiza- 
tion problem. The discussion that follows is somewhat different, and is derived from the refereed 
quantum games framework of [GW07J. 

Suppose hereafter in this section that <& = Oo — Oi for quantum channels Oo, Oi 6 T [X, y), 
and consider the semidefinite program whose primal and dual problems are as follows: 



Primal problem 



maximize: 
subject to: 



<J(0),W) 
W < ty ® p, 

w g Pos (y (g> X) 



Dual problem 
minimize: ||Try(Z) H^, 
subject to: Z > J(<J>), 

z g Pos (y - 



) X) 



As in the previous section, these problems can be matched to the formal description of a semidef- 
inite program (Y, C, D), for which strong duality is easily proved. Our goal will be to prove that 
the optimal value of this semidefinite program is given by \ ||| <3E> ||| ^ - 

Given that O is the difference between completely positive super-operators, it holds |GLN05 , 

IRW05llWat05ll that 



IOI 



max 

ueS(x®x) 



(*®1 LW )(UM 
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Given that the operator (O ® lux) ) ( uu * ) is the difference between two density operators for every 
u G S{X ® W), it follows that 

IIIOIHj = 2max{^P,(0®l L((Y) )(uu*)^ : u G «S(# <8>#), P G Pos (y g) AT) , P < ly®x] . 

Now, for every unit vector u G X ® X there is a corresponding operator B G L (X) with 
|| B || 2 = 1 such that 

u= (Ei,j,B)ei®ej. 

l<i,/'<« 

For this choice of £> we have 

(ly ® B)/(*)(ly ® B*) = (O ® 1 LW )(««*) 

and so 

lllOlll! = 2 max ((1 ® B*)P(1 ® B)J(O)) 

where the maximum is over all B G L (X) with || B || 2 = 1 and P G Pos (3^ ® A") with P < ly®^- 
Now define sets Q and 1Z as follows: 

Q = {R £ Pos (y ® A") : R < ly <8> p for some pGD(Af)} ( 

7i = {(ly®B*)P(l 3 ;®B) : B G L(X), P G Pos (;V ® AT), ||B|| 2 = 1, P<ly®x}- 
It holds that 

Hl^lll! =2 sup <J(<&), X) 

Xell 

while the optimal value of the semidefinite program is 

« = sup </(*), X). 

The fact that a. = j therefore follows from the equality Q = 1Z, which is easily proved by 

selecting p or B so that p = B*B. 

5 Connections with known results 

This section describes two interesting connections between the semidefinite programming formu- 
lation from Section|3]and known results, the first being directly about completely bounded norms, 
and the second concerning the fidelity function. 

Spectral norms of Stinespring representations 

The following theorem gives an alternate characterization of the completely bounded trace norm 
(or diamond norm). Proofs can be found in Kitaev, Shen and Vyalyi [KSV02J and Paulsen |Pau02J. 
The two proofs use rather different techniques, and here the theorem is proved in a third way 
using semidefinite programming duality. 

Theorem 5. For every super-operator O G T {X , y), it holds that 

|||0||| 1= inf ||AL||BL, (4) 

(A,B) 

where the infimum is over all Stinespring pairs (A, B)for <E>. 
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Proof. For any Stinespring pair (A, B) of O, where A,B 6 L {X,y ® Z), and for any choice of 
W = C fc , it holds that 



(<D®1 L(VV) )(X) 



= ||Tr z [(A(g)]L w )X(B*(g)]L vv )]|| 1 
< || (A <g> l>v)X(B* <g> l>v) \\ x 



< IIAI 



l x lli ll B l 



for all X E L(Af ®W). It follows that IIIOH^ < pl^. 

To prove that the infimum is no larger than 1 1| <3> 1 1| 1 , first choose an arbitrary Stinespring pair 
(A, B) of O, where A,B 6 L [X , y <8> Z). The optimal value for the dual problem stated in Sec- 
tion|3]does not change if Z is restricted to be positive definite, provided we accept that an optimal 
solution may not be achieved. We therefore have 

\\\<S>\\\ 2 1 =mf{\\A*(ly®Z)A\\ !X> : ty®Z>BB*, ZePd(Z)}. 

Thus, for a given e > 0, we may choose Z E Pd (Z) such that 

||A*(l y ®Z)A|L<(|||0||| 1 + £ ) 2 

and ty <S> Z > BB* . This second inequality is equivalent to 



So now we have that 



and it holds that 



y 



'1/2 



Z -l/2\ BB * ( 1 



yy 



-1/2 



A 



r-l/2 



< 1. 



< lll^llll+ £ ' 



ty®Z l/1 ) A, [ly^Z' 1 - ) B 



r-l/2 



is a Stinespring pair for O. This establishes that the infimum equals IH^HIi in the expression HJ, 
which completes the proof. □ 

Connection with fidelity 

Consider the semidefinite program from Section [3l for the special case where X = C. Replacing 
A and B with vectors u,v $L y ® Z , and making simplifications, the problems become as follows: 



Primal problem 
maximize: (vv*,W) 
subject to: Trj;(W) < Try (uu*) , 

w e Pos iy®z). 



Dual problem 
minimize: (Try(uu*),Z) 
subject to: ty <g) Z > vv* , 
Z E Pos(Z). 



The quantity that is represented by the optimal value of these problems is given by the fidelity 
function, which is defined as 



F(P,Q) 



vWQ i = Tr V VpQ^P 
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for positive semidefinite operators P and Q. In particular, the optimal value (for the primal and 
dual problems) is 

F(TTy(uu*),Tr y (vv*)) z , (5) 

as is now explained. 

First, the optimal value of the primal problem follows from Uhlmann's Theorem |Uhl76| , 
which is as follows. 

Theorem 6 (Uhlmann's Theorem). Let y and Z be finite-dimensional complex vector spaces, and let 
P, Q G Pos (Z) be positive semidefinite operators, both having rank at most dim(3^). Then for any choice 
ofvey^Z satisfying Try(vv*) = Q, it holds that 

F(P,Q) =max{|(u,y)| : u G y ® Z, Tr y (uu*) = P} . 

It is straightforward to obtain from this theorem (along with simple properties of the fidelity) the 
following corollary, which is precisely the statement that the optimal primal value of our semidef- 
inite program is given by the fidelity. 

Corollary 7. Assume u,v G y (g) Z are vectors, and let P = Try(uu*) and Q = Try(vv*). Then 

F(P,Q) 2 = max {(vv*,W) : W G Pos (y ® Z) , Tr y (W) < P} . 

The optimal value of the dual problem is, of course, equal to © by strong duality. A different 
way to evaluate the optimal dual value begins with the following simple proposition. 

Proposition 8. For any vector v G y <8> Z and any positive definite operator Z G Pd(Z) it holds that 
l y ®Z> vv* ifandonlyif{Tr y (vv*),Z- 1 ) < 1. 

Proof. It holds that ly <g) Z > vv* if and only if 

% ® Z~ 1/2 ) vv* (ly ® Z~ 1/2 ) < l y ® z . (6) 



Given that the operator on the left-hand-side of (© is positive semidefinite and has rank equal to 1, 
we have that © is equivalent to 



ty®Z- 1/2 



< 1, 



-1/2 



< 1. 



-1 



which in turn is equivalent to 

Tr (jl y ®Z- 1/2 ^ vv* (ly§ 

As 

Tr ((ly ® Z~ 1/2 ) vv* (t y ®Z- 1/2 ^ = (lr y {vv*),Z 

the proof is complete. □ 

We have that the optimal value of the dual problem does not change if Z is optimized over only 
positive definite rather than positive semidefinite operators (again accepting that the optimal 
value may not be achieved for such an operator). Combined with the proposition just proved, 
we find that the optimal dual value is given by 

j S = inf{(Try(ww*) / Z) : ZG Pd(Z), {Tryivv*)^- 1 ) < l} . 

That this value is given by (O follows from a different characterization of the fidelity due to Alberti 
HAlb83L 
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Theorem 9 (Alberti's Theorem). Let P,QG Pos (Z) be positive semidefinite operators. Then 

(F(P,Q)) 2 = inf (P / Z)(Q / Z- 1 ). 

We have therefore established a simple and precise sense in which Uhlmann's Theorem and Al- 
berti's Theorem are dual statements in finite dimensions, each implying the other. 
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